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Abstract: Obtaining ultracold samples of dipolar molecules is a current challenge which 
requires an accurate knowledge of their electronic properties to guide the ongoing 
experiments. Alkali hydride molecules have permanent dipole significantly larger than 
those of mixed alkali species and, as pointed out by Taylor- Juarros et at [Eur. Phys. J. D 31, 
213 (2004)] and by Juarros et al. [Phys. Rev. A 73, 041403 (2006)], are thus good candidates 
for molecule formation. In this paper, using a standard quantum chemistry approach based 
on pseudopotentials for atomic core representation, large Gaussian basis sets, and effective 
core polarization potential, we systematically investigate the electronic properties of the alkali 
hydrides LiH to CsH, in order to discuss general trends of their behavior. We computed (for 
the first time for NaH, KH, RbH, and CsH) the variation of their static polarizability with 
the internuclear distance. Moreover, in addition to potential curves, we determine accurate 
values of permanent and transition dipole moments for ground and excited states depending 
on the internuclear distance. The behavior of electronic properties of all alkali hydrides is 
compared to each other, in the light of the numerous other data available in the literature. 
Finally, the influence of the quality of the representation of the hydrogen electronic affinity in 
the approach on the results is discussed. 

Key words: LiH, NaH, KH, RbH, CsH, alkali hydrides, permanent dipole and transition 
dipole moments, static dipole polarizabilities. 

PACS Nos.: 31.15.ac,31.15.ap,31.50.Bc,31.50.Df 

1. Introduction 

Alkali hydrides have continuously attracted the interest of researchers in various areas. Indeed, they 
are among the molecules with the simplest electronic structure, allowing very detailed comparisons 
between different theoretical models, and with available experimental results. The astrophysical and 
cosmological relevance of LiH is well known [1—3], as it is believed to be formed in the early universe 
by radiative association of H and Li [1]. The determination of the abundances of chemical elements in 
stellar atmospheres and the modeling of the chemistry of these elements in the early universe requires 
a knowledge of the cross sections for various collisional processes such as inelastic collisions between 
hydrogen and alkali atoms (Li, Na) [4,5], radiative charge transfer between alkali atoms and H + [6], 
dissociative recombination [7], or neutralization collisions between alkali ions and H~ [8,9]. Molecular 
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reaction dynamics of alkali atoms and molecular hydrogen involve alkali hydrides as products and a 
good description of these systems is required to analyze the reaction mechanism [10, 11]. Cote et 
al. [12, 13] suggested that the formation of Bose-Einstein condensates of a dilute gas of hydrogen 
could be achieved by cooling hydrogen atoms in a buffer gas of alkali atoms, which critically depends 
on the related scattering lengths. Recently Juarros et al. have predicted that formation of ultracold LiH 
or NaH molecules is possible using various processes such as stimulated Raman photoassociation or 
stimulated one-photon photoassociation [14-16]. Due to its large permanent dipole moment, the LiH 
molecule [17], and possibly the NaH molecule, may be considered as candidates for Stark deceleration. 
Also, following the ideas of ref. [18] related to alignment and orientation of polar molecules under 
the influence of combined static and laser fields, alkali hydrides could be good candidates for such 
experiments. An accurate knowledge of the electronic properties of these molecules are prerequisite to 
the extension of these investigations to heavier alkali hydrides. 

Following previous work on alkali dimers in our group, we performed systematic computations of 
the electronic properties of alkali hydrides from LiH to CsH (Section 2) as well as the monocations 
formed from these molecules. New results for the dependence on internuclear distance and on the vi- 
brational index of their dipole polarizabilities are presented in Section 3. We also compare our results 
with numerous previous theoretical data for potential curves and permanent and transition dipole mo- 
ments, identifying general trends for these quantities (Section 4), and allowing for an analysis of their 
accuracy. In particular, the issue of the inclusion of a correction for the hydrogen electron affinity is dis- 
cussed. Most of the computed results of the present work are available as supplementary material pro- 
vided by the editor. In the following we will use atomic units for distances (a =0.0529 17720859 nm) 
and energies (2R oc = 219474.63137 cm -1 ) except where otherwise stated. 

2. Computational details and potential curves for neutral and ionic systems 

As in our previous studies on alkali dimers, the alkali atoms are described by ^-dependent pseu- 
dopotentials for the ionic cores [19, 20] including effective core polarization terms [21, 22], and by a 
large set of unconnected Gaussian functions for the valence electron. For all atoms but Cs we used 
the Gaussian basis sets labeled as "B" in ref. [23,24], while for Cs we use the basis set labeled as 
"B"' in ref. [25]. We also use the same core-polarization potentials (CPP) than in our previous work, 
with the reported values for the cut-off radii. For the H atom we set up a new large [10s5p2d] basis 
set of uncontracted Gaussian functions with the following exponents: 195.5, 27.6, 6.3, 1.8, 0.5983, 
0.22, 0.09, 0.04, 0.019, 0.006 for s basis functions; 1.05, 0.22, 0.0669, 0.024270, 0.00669 for p ba- 
sis functions; 0.34, 0.061 for d basis functions. This choice globally improves the calculated atomic 
level energies compared to previous works. We obtain the following discrepancies: AEi s — 3.2 cm -1 , 
AE 2s = 2.1 cm -1 , AE 2p = 7.4 cm -1 , AE 3s = 29 cm -1 , AE 3p = 51 cm -1 . For instance, Huzinaga 
et al. [26] used a large basis set of contracted Gaussian orbitals which gave exact energies for the Is, 
2s, 2p, 3s and 3p levels. In their recent series of papers, Gadea and coworkers [27-30] designed a 
[7s3p2d] Gaussian basis set contracted to [5s3p2<i], close to the one of Geum et al. [31] who used a 
[7s3p2rf] contracted to [6s3p2d] basis set. The former authors obtained deviations of A_Ei s ^11 cm -1 
and AE 2s ~ 700 cm" 1 [27], and the latter AE ls ~ 64cm -1 and AE 2s ~ 35 cm -1 . Moreover our 
basis also provides an improved hydrogen electronic affinity with a deviation of about 140 cm -1 from 
the exact value of ref. [32], i.e. three times smaller than the value quoted by Boutalib and Gadea [27] 
but slightly larger than the value of Geum et al. 111 cm -1 ). 

The molecular calculations are performed just like in our previous papers [23,25, 33]. Molecular 
orbitals are determined by a restricted Hartree-Fock calculation including core polarization, yielding 
electronic energies and wave functions of the molecular ion. Electronic energies and wave functions of 
the neutral system are deduced from a full configuration interaction (CI) of two-electron configurations, 
through the CIPSI quantum chemistry code of the quantum chemistry group in Toulouse (France) [34]. 
Our procedure automatically provides potential curves for the related alkali hydride ions treated as 
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effective one-electron systems, and it is worthwhile to examine them as a first check of its accuracy 
compared to other approaches. It is striking to see in Table 1 that even for the ground state of such 
simple systems, published values vary significantly for the binding energy of LiH + and NaH + , while 
only few data are available up to now for the heaviest species KH + , RbH + , and CsH + . To our knowl- 
edge only one experimental value has been published up to now, namely the CsH + binding energy 
extracted from differential cross section measurements [35]. Very few measurements have been re- 
ported for other molecular states. We also list in Table 1 data for the first excited state A 2 H + which has 
been widely studied in the context of the theoretical calculation of cross sections for charge exchange 
reactions [1,6, 36-39] . 

Most of the calculations of alkali hydride ion potential curves have been carried out with effective 
core potentials. In a series of recent papers, Magnier [4CM-2] extended the work of Alikacem and 
Aubert-Frecon [43] by including effective core polarization potentials to their Klapisch model potential 
approach [44]. Various kinds of pseudopotentials have been used by several authors [36,37,39,45-48], 
the closest approach to ours being the work by Berriche and Gadea [49] for LiH + . A few ab initio 
all-electron calculations have also been performed many years ago, by Rosmus and Meyer [50] with 
the coupled electron pair approach, by Vojtik [51] using the MRDCI approach, and by Dalgarno et 
al. [1] and Olson et al. [52] with SCF and CI approaches, for the treatment of radiative association in 
an astrophysical environment. It is tedious to achieve an enlightening comparison of our data with these 
works, as most of them are quite old and were probably limited by computational facilities at that time. 
Moreover, several calculations reported binding energies D e for the A excited state only, with which 
our results reasonably agree, but without reporting similar data for the ground state. Our values for the 
LiH + ground state agree well with those of Magnier [40], but are slightly larger than those obtained 
by the same author on NaH + and KH + [41,42]. For the A state the present D e values for LiH + and 
NaH + are almost identical to those predicted by Magnier [40,41] but slightly smaller for KH + [42]. 
Also, the depth of the potential well of the A state of CsH + lies within the error of the data reported in 
ref. [35]. 

Next, the accuracy of our full CI results yielded by our treatment of the neutral molecule as an 
effective two-electron system, can be analyzed from Tables 2 and 3, which display the equilibrium 
distances R e , the harmonic constant u e and the potential well depth of the alkali hydrides electronic 
states together with available experimental values and with several other recent theoretical results. We 
concentrate our analysis on the fourth lowest 1 S+ and two lowest 1 II electronic states, as they can be 
reached by electric dipole transition from the ground state. However, we also performed calculations 
for triplet states, and all these results are collected in the supplementary material attached to this paper. 
We label the molecular states according to the standard spectroscopic notation, i.e. X, A, C, E corre- 
sponding to the four lowest 1 X + states, and B = (l) 1 !!. Similarly, triplet states are labeled with the 
lowercase indexes a, b, and c, for the (1) 3 S + , (1) 3 II, and (2) 3 £+, respectively. In refs. [27-30], the 
authors labeled the (4) 1 S + as D state. Our convention for the sign of the permanent dipole is that a 
positive value corresponds to the charge distribution A + H , where A is an alkali atom. 

The extensive review by Stwalley et al. [54] addresses most of the works already published on 
alkali hydrides. In particular, this review reports on experimental studies of the first 1 £+ state (the A 
state) from NaH to CsH. Vidal et al. [55,56] investigated the spectroscopy of A and B (the lowest 1 II) 
states, and Yang et al. [57] and Hsieh et al. [58] the A state in KH and in CsH, respectively. 

We briefly recall below some of the theoretical calculations, mainly those yielding tables of per- 
manent or transition dipole moments, for further reference in the next sections. Ab initio multicon- 
figuration self-consistent-field calculations with configuration interaction (MCSCF-CI) of permanent 
and transition dipole moments of LiH have been reported by Docken and Hinze [59] using a Slater- 
type basis set, in contrast with Gaussian basis sets employed in the present calculations. Partridge 
and Langhoff [60,61] repeated such computations with an extended Slater-type basis set. Roos and 
Sadlej [62] used the complete active space SCF (CASSCF) approach to calculate the permanent dipole 
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moment and the polarizabilities of the ground state as functions of the interatomic distance. Sachs et 
al. determined the electronic properties of several states of NaH [63,64]. Langhoff et al. extended 
their previous investigations to NaH, KH, and RbH [65] using near Hartree-Fock quality Slater basis 
sets and incorporating electron correlation through a coupled-pair formalism. Laskowski et al. [66] 
studied the two lowest 1 S+ states in CsH, using a similar approach to ours: effective ^-dependent 
pseudopotentials derived from relativistic Hartree-Fock calculations and including core polarization 
terms for atomic core representation, and Gaussian basis sets in the CI calculation. Carnell et al. [67] 
investigated the first seventeen states of CsH with ab initio multi-reference configuration interaction 
(MRDCI) calculations. Combining existing experimental data and ab-initio calculations, Zemke et al. 
determined the dipole moment for the A 1 ^ -> B l Ii -> and B l Ii -> tran- 

sitions and related radiative transition probabilities in LiH [68, 69], and the permanent dipole moment 
for the X and A states in NaH [70]. Finally, Camacho et al. [71] built the -> transition 

dipole moment of KH from available experimental data. 

More recently Gadea and coworkers have performed a series of studies of the adiabatic and diabatic 
molecular states of alkali hydrides [27-30]. As in the present work, the authors used the ^-dependent 
pseudopotentials of Durand and Barthelat [19,20], and core-valence correlation from ref. [21,22]. Fur- 
thermore, an i?-dependent correction of the ion-pair diabatic curves is introduced to account for basis 
set limitations and to ensure an improved value for the ground state energy. Permanent dipole mo- 
ments for several 1,3 S + states and transition dipole moments between 1 £ + states have been displayed 
for KH [28], RbH [29], and CsH [30]. Khelifi et al. have determined the radiative and nonradiative 
lifetimes of the and C 1 S+ vibrational levels of the KH molecule [72,73]. 

In Table 2, we see that theoretical values for R e are systematically slightly smaller than the exper- 
imental ones, especially for CsH. This could be due to the absence of short-range repulsion terms in 
the ECP calculations. Our resulting values for D e are slightly larger than similar calculations of Gadea 
and coworkers (column (a) in the Table) [27-30], which is a manifestation of the influence of the large 
basis we used in our work. However it is well known that ECP-type calculation is not a perfect vari- 
ational approach, so that it may happen that the computed well-depth exceeds the experimental one, 
especially when large ionic cores are involved (see the CsH case). This pattern is also visible in results 
from Dolg [74]. They rely on fully relativistic pseudopotentials, and molecular calculations are based 
on Dirac-Hartree-Fock and CI approaches. Core-polarization terms have been introduced, as well as 
corrections for core-core repulsion. Equilibrium distances are also found shorter than the experimental 
ones, while D e are larger than the experimental ones for all species. When an empirical correction 
for the hydrogen electronic affinity is included by Gadea and coworkers for KH, RbH, and CsH (see 
section 4), the values for D e are increased and become closer to the experimental ones. However, such 
a correction cannot be considered as a variational procedure, so that their D e value for CsH exceeds 
the experimental one by about 462 cm -1 . We will comment more about the effect of their ion-pair 
correction on dipole moments in section 4. 

We summarized the main spectroscopic constants of the four lowest 1 S + excited states, and of the 
lowest : n state in Table 3, compared to experimental observations, and to selected theoretical predic- 
tions. As already pointed out, Gadea and coworkers [27-30] have investigated all systems but NaH and 
their papers include extensive comparisons with other predictions which are not duplicated here. For 
NaH our values are compared to those of Lee et al [76], who used the large-core pseudopotential of the 
Stuttgart group [77] and accounted for core polarization following ref. [21]. SCF and CI calculations 
have been performed with the MOLPRO package. The X and A states of KH and RbH have been 
investigated by Garcia et al [80] using a small-core pseudopotential [77], and the difference-dedicated 
configuration interaction (DDCI) method to account for core- valence correlation effects. Our values 
are generally in good agreement with those of Gadea and coworkers although the Gaussian basis sets 
used in those work for describing the alkali and hydrogen atoms differ from ours. It is difficult to know 
if the somewhat large differences with the values for NaH of [76] are related to the differences in the 
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pseudopotential or to those in the basis sets. 

3. Static dipole polarizabilities 

If the z axis is chosen along the internuclear axis in the molecule-fixed reference frame (x,y,z), 
there are two independent components of the molecular polarizability tensor, i.e., the parallel compo- 
nent a|| = a zz and the perpendicular one a± = a xx = a yy . Two related quantities are usually defined: 
the average polarizability a — (a\\ + 2a±)/3 and the polarizability anisotropy 7 = an — a±. As in 
our previous work on the static dipole polarizabilities of alkali dimers [82], we determined the static 
dipole polarizabilities of alkali hydrides in their ground state with the finite-field method [83], using 
electric fields between 0.0003 and 0.0005 a.u. to remain in the perturbative regime. We calculated both 
components of the static dipole polarizability for all systems as functions of the internuclear distance 
R (Figure 1). For all systems, the i?-variations of an and a± are similar, with an increase of ay with 
increasing alkali mass, and a smooth increase of a± with R. To our knowledge, no other i?-variation of 
the alkali hydride dipole polarizability has been previously published, except for LiH in ref. [62] which 
displays a dependency in reasonable agreement with the present one (Figure 1(c)). Looking at large R, 
we note however that the lithium polarizability seems to be slightly underestimated in ref. [62]. Also, 
Kolos and Wolniewicz [84] found similar variations of the H 2 polarizabilities than the present ones, 
but with smaller values, as expected. 

By integrating the polarizabilities over the vibrational wave functions, we deduce the u-dependency 
of a and of 7, which regularly increase from LiH to CsH, just like for the series of alkali pairs (Figure 
2(a)). The anisotropics reach quite large values (up to 967 au. for CsH at v = 24, i.e. two times 
larger than the cesium atomic polarizability), magnifying the influence of the hydrogen atom when the 
molecule is vibrationally excited, in contrast with the v — case. No other values have been published 
in the past. We also display in Figure 1(b) the permanent dipole moment of the alkali hydride ground 
state which show similar w-dependence for all the systems, with a maximum followed by a rapid 
decrease to zero, as expected from comparable calculations for alkali pairs [23]. 

In our previous work about the ground state polarizabilities of all homonuclear and heteronuclear 
alkali diatomics [82], both an and a± components were found to vary linearly with an effective volume 
V e ff = 47r(i? e ) 3 /3, a || varying two times faster than a± with V e ff. This suggested the picture of an 
effective elliptic charge distribution for the alkali pair at the equilibrium distance. In alkali hydrides, 
this scaling law remains valid only for an, while a± does not monotonically vary along the series of 
alkali hydrides. This may be due to the weakness of the hydrogen polarizability, which is at least one 
order of magnitude smaller (4.5 a.u. [85]) than the alkali ones ( [82]). Therefore the effective elliptic 
charge distribution of the molecule is progressively dominated by the polarizability of the alkali atom, 
along the series from Li to Cs. 

Most previous theoretical works reported polarizability values at the experimental equilibrium dis- 
tance Rl xp of the ground state, which are assembled in Table 4, together with the present results for 
a II , a±, a and 7 calculated either at R% xp [54] or for v — 0. The LiH molecule is the most extensively 
studied system, and we give only the most significant values. Additional references may be found in 
the quoted papers. Coupled Hartree-Fock perturbation theory has been used by Lazzeretti et al. [86], 
and multiconfigurational time-dependent Hartree-Fock (MCTDHF) method by Sasagane et al. [87]. 
Multiconfiguration self-consistent field (SCF) approaches have been used at various levels by several 
authors [62,88-90]. In these papers, the finite field method [83] or the charge perturbative approach [89] 
have been used to extract polarizabilities. A time-dependent gauge invariant method (TDGI) has been 
used by Rerat et al. [91], while Vrbik et al. [92] employed a diffusion quantum Monte-Carlo (DQMC) 
approach. Sadlej and coworkers are the only authors who investigated static dipole polarizabilities for 
other species than LiH using many-body perturbation theory (MBPT) and coupled-cluster theories 
combined with finite field methods [93,94]. Most of the above-quoted authors presented predictions 
obtained with different basis sets or at different levels of approximation, and we only report in Table 4 
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Fig. 1. (a) a || and (b) a± components of the static dipole polarizability of the alkali hydride ground state as 
functions of the internuclear distance, (c) For LiH, our results for ay (full line) and a± (dashed lines) are compared 
with those of Roos and Sadlej [62] (X-es). 
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values obtained with their most sophisticated model. 

Values for the polarizabilities at R e are similar to those for v = for all species but CsH within 
2%, simply because the steepness of the polarizability functions around R e is compensated by the sym- 
metric character of the v = wavefunction. For CsH, the 9% difference comes from our equilibrium 
distance, which is slightly smaller than the experimental one. For LiH and NaH our predictions agree 
well with other values with the exception of the values of Gready et al. [88] obtained with a too small 
basis set. The relatively wide variation of the theoretical values clearly emphasizes the high sensitivity 
of the polarizability to the size of the basis sets and to correlation effects. Our predictions of ay for 
heavier systems lie between the two quoted values of ref. [93], while our values for a± are slightly 
smaller. Note that the dipole moments calculated with their approximation (b) in that paper are closer 
to our values than their (a) values. 

We are not aware of any polarizability measurement in alkali hydrides. An old experimental value 
of the polarizability anisotropy for LiH with a large error bar (7 = 1.7 ± 4. au) has been quoted by 
Stevens and Lipscomb [95] as an unpublished result from Klemperer et al. 

4. Trends of the permanent and transition dipole moments of alkali hydrides 

We report in Table 4 our values for the permanent dipole moment [i z at the experimental equilib- 
rium distance of the ground state. We note that they are quite a bit larger than those of alkali dimers [23] 
(ranging between 0.5 and 5 Debye), due to the large asymmetry of the charge distribution. Our val- 
ues show good agreement with the limited experimental data available (currently only for LiH and 
NaH [96, 97]). Table 4 reflects the situation which held until recently for the theoretical works: many 
of them were devoted to LiH, and only a few to heavier species. In particular, we obtain a very good 
agreement with the most elaborate model of Sadlej and Urban [93]. 

We also computed the i?-variation of the permanent and transition dipole moments for a selection 
of electronic states. It is well-known that such functions represent sensitive tests for the accuracy of the 
computed electronic wave functions. As mentioned earlier, many theoretical works have been already 
published for these quantities, and we concentrate on systematic trends in their behavior along the 
series of alkali hydrides. It is worth noting that such an analysis can be performed due to the wealth 
of data concerning alkali hydrides, in contrast with the few data available for alkali dimers (see for 
instance ref. [33]). 

In the following figures, we display the computed i?-variation of the permanent dipole moments 
for the X, A, B, C, a, and b states of alkali hydrides, compared to other available studies. As expected 
the most studied state is the X ground state (Fig. 3(a)). Its permanent dipole moment shows similar 
regular variations with R for all species, reaching a higher maximum amplitude at larger distance with 
increasing mass and polarizability. The agreement between all methods is generally satisfactory, with 
the exception of ref. [67] where the maximum value is found to be about 20% larger. The permanent 
dipole moments for the A state (Fig. 3(b)) vary like the one for the X state, reaching their maximum 
value at larger distances. Several other determinations exist in the literature for the LiH and NaH A 
state. We found a remarkable agreement for LiH with refs. [59,60]. In contrast, the maximum value 
found by Sachs et al. [63] is almost two times smaller than ours. As for the ground state, the maximum 
value found by Carnell et al. for CsH is larger (by about 25%) than ours, while the one by Laskowski et 
al [66] is about 25% smaller than ours. To our knowledge, the only other calculations for KH, RbH, and 
CsH have been reported by Gadea and coworkers in refs [28], [29], and [30], respectively. No tables of 
numerical data are available in these papers, but a satisfactory agreement is visible, as long as the eye 
can judge from their figures (see the discussion for RbH later in this section). Finally our calculations 
for both the X and A states of NaH are in excellent agreement with those of Zemke et al. [70] who 
extracted the dipole moment functions from an analysis combining ab-initio and experimental results. 

As expected, the i?-variation of the A — X transition dipole moment (Figure 4a) reflects both varia- 
tions of the X and A permanent dipole moments through the corresponding electronic wave functions: 
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two extremum are visible, whose positions correspond to those obtained for the X and A permanent 
dipole moment functions. The same figure also displays results concerning the transitions B — X and 
B — A involving the lowest B 1 !! excited state (Figure 4b, c). An excellent agreement is found again 
with refs. [59, 60] for all these transitions in LiH. Quite unexpectedly, the agreement with the values 
of Sachs et al [63] for NaH looks better than for the A permanent dipole moment. For CsH, the curve 
by Laskowski et al [66] is close to ours in the region of the maximum value of the X dipole mo- 
ment (below 10 a.u.), and not in the region of the one of the A dipole moment (around 15 a.u.). As 
in the previous figure, results from ref. [67] disagree with ours. Again, the agreement for the A — X 
transition dipole moments seem satisfactory when looking at figures from refs. [28-30] for KH, RbH, 
and CsH. While not displayed here, results are also available for transition dipole moments among the 
lowest triplet X+ and II states (c — a, b — a, c — b) in LiH [59] and NaH [63], in good agreement 
with the present results. Experimental values are only available for A - X and B — X transitions in 
LiH [68,69], and our values are in excellent agreement with them. A determination of the A — X tran- 
sition dipole moment variation in KH has been recently extracted from the experiment by Camacho et 
al. [71] (which is to our knowledge the only experimental data for this molecule) showing a reasonable 
agreement with our value. 

In Figure 5 we display the permanent electric dipole moment for higher excited singlet states (B 
and C), and for the two lowest triplet a and b states. While less data are available for comparison, 
we observe the same trend than in the previous figures: the systematic investigation of Carnell et al. 
ref. [67] really seems to miss most of the character of the excited electronic wave functions of CsH, 
while our results are in agreement with those of ref. [60] for the B state of LiH, and of ref. [63] for the 
B, a, and b states in NaH. A visual satisfactory agreement is also observed with the figures displayed 
in refs. [28-30]. 

As expected, all of these dipole moments vanish with increasing distance, with a smaller rate with 
increasing alkali mass. The C state is the most excited state displayed in the figure - correlated to the 
second excited ls+3s, ls+4s, ls+5s, ls+Ad, and ls+5d dissociation limits for LiH, NaH, KH, RbH, 
and CsH, respectively- so its magnitude reaches larger values than for the other states. For all alkali 
hydrides, the C dipole moment exhibits irregular variations, which is the well-known manifestation of 
the contribution of the ion-pair configuration in the electronic wave function (see Figure 6): as long 
as this ion-pair character dominates the wave function (between the avoided crossings marked with 
circles in Figure 6), the permanent dipole moment can take large values, until the distance where the 
covalent character shows up again. In LiH and NaH, the presence of the avoided crossing C and E 
states at short distances is also visible in Figure 5c (see arrows). Quite unexpectedly, the CsH dipole 
moment is not the largest one of the series, as the ion-pair character cannot develop towards distances 
larger than 22 a.u. due to the somewhat smaller energy spacing between the excited dissociation limits, 
compared to the other species. 

The influence of the ion-pair character is also visible in the transition dipole moments among these 
excited electronic states, as visible on Figure 7 for the C — A and E — C transitions. The transition 
dipole moment functions change very abruptly with the internuclear distance (see arrows), which is 
the expected manifestation of local change of the electronic wave functions at the location of avoided 
crossings between the related potential curves. 

The treatment of the ion-pair contribution is actually a sensitive issue for discussing the quality 
of the representation of the electronic wave functions of the molecules in various approaches. Indeed, 
such a character is not explicitly introduced in the standard Gaussian basis sets for neutral atoms, and 
would need in principle an infinite number of basis functions to be accurately represented. In the present 
work, we set up a large basis for the hydrogen atom, in order to reproduce the atomic energy levels, and 
hopefully to take in account as much as possible of the ion-pair character. An indication for this is our 
computed value of the hydrogen electronic affinity EA = 5938 cm -1 , compared to the experimental 
value EA exp = 6083 cm -1 [32]. Let us note that in ref. [76] mentioned in Table 3 the authors used a 
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Fig. 4. Computed transition dipole moments for (a) the yl 1 E + — X X E + transition, (b) the B 1 !! - X X E + 
transition, (c) the B 1 !! — A 1 E + transition, as functions of the intemuclear distance. Our results are: LiH (black 
full line), NaH (red dashed line),KH (green dot-dashed line), RbH (blue dotted line), CsH (orange long-dashed 
line). Other determinations are displayed from ref. [59] (black full circles), ref. [60] (black open circles), ref. [63] 
(open squares), ref. [68,69] (plus signs), ref. [71] (crosses), ref. [66] (full triangles), and ref. [67] without or with 
relativistic effects (dotted line and full circles, or dotted line with open circles, respectively). 
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Fig. 5. Computed permanent dipole moments of B 1 !!, C 1 E + , a 3 E + , and fe 3 II states as functions of the internu- 
clear distance. Our results are: Our results are: LiH (black full line), NaH (red dashed line),KH (green dot-dashed 
line), RbH (blue dotted line), CsH (orange long-dashed line). Other determinations are extracted from ref. [59] 
(black full circles), ref. [60] (black open circles), ref. [63] (open squares), and ref. [67] without relativistic effects 
(dotted line and full circles). Arrows for each species correspond to the location of avoided crossings in the related 
potential curves in Figure 6. 
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Fig. 6. Potential curves of the alkali hydrides E + states in the region where the ion-pair contribution dominates 
the electronic wave function. The A, C, and E states are drawn with dot-dashed, full, and dashed lines respectively. 
The origin of energies is taken at the double ionization limit of the atoms at infinite distance. Avoided crossings 
marked with circles correspond to the features discussed in Figure 5c and in Figure 7. 
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Fig. 7. Computed transition dipole moments for (a) the C^E 4 " - A X E + transition, (b) the _E 1 E + - C^E 4 " 
transition as functions of the internuclear distance. Our results are: Our results are: LiH (black full line), NaH (red 
dashed line),KH (green dot-dashed line), RbH (blue dotted line), CsH (orange long-dashed line). We reported in 
panel (a) the value from ref. [67] without relativistic effects (dotted line and full circles). Arrows for each species 
correspond to the location of avoided crossings in the related potential curves in Figure 6. 
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7s4p2d [6s4p2d] Gaussian basis set [98] for the H atom, yielding a good determination of the hydrogen 
affinity with deviation from the experimental [32] of ~ 145cm -1 ,Le. similar to our value, while the Is 
and 2s states are less accurate than in our work, with A ~ 64 cm -1 and AE <~ 77 cm -1 respectively. 

The value of EA determines the location of avoided crossings between potential curves with cova- 
lent and ion-pair character. In a series of papers, Khelifi et al. [28,29] and Zrafi et al. [30] addressed this 
issue in detail in KH, RbH, and CsH, respectively. We note that the overall agreement with the numeri- 
cal values for RbH from Khelifi [99] is excellent (Figure 8), which is not surprising as our methods are 
quite close. Using a diabatization procedure, they were able to build an effective Hamiltonian which 
separates the covalent and the ion-pair characters of the electronic states. Then they account for their 
discrepancy on the hydrogen EA (405 cm -1 in their case) as an empirical negative correction added 
to the corresponding diagonal element of the effective Hamiltonian. The related avoided crossings in 
the adiabatic potential curves are then shifted towards larger internuclear distances, as well as the as- 
sociated abrupt variations in the transition dipole moment functions (see dashed lines in Figure 8 for 
the case of RbH). In contrast, we see from our calculations that the increase of the basis set size, which 
yields a value for the hydrogen EA as close as 145 cm -1 from the exact value, only provides a small 
shift of these patterns toward large distances compared to the results of ref. [29] without the empirical 
correction. As tentative interpretation, we estimate that the position of the avoided crossing between 
the E and the C curves around 26.5a (Figure 8c) would be shifted only by about 0.4ao if we would 
lower our E potential curve down by 145 cm -1 . However in our results, the associated peak in the 
C — E transition dipole moment around 26.5a (Figure 8c) is separated by about 0.8a from the one of 
ref. [29] when the empirical correction is included. Therefore, it is not obvious if such a global energy 
shift of the ionic configuration accurately describes the ionic configuration as this correction does not 
change the electronic wave function of the negative ion, as also illustrated for instance by the excessive 
increase of the CsH potential well depth by 462 cm -1 . 

5. Conclusion 

In this paper, we computed the electronic properties of alkali hydrides from LiH to CsH, includ- 
ing potential curves, permanent and transition dipole moments, and static dipole polarizabilities. We 
performed a systematic investigation which allowed us to ensure a similar numerical accuracy for all 
species, in the same spirit as the numerous studies published by Gadea and coworkers. Our results are 
found in good agreement with available experimental data, so that our systematic computations for 
all species are useful to estimate the accuracy of other available theoretical results. Apart for LiH, we 
determined the variation of the polarizabilities with internuclear distance for the first time, and their 
values for vibrational levels still await for experimental confirmation. Furthermore, general trends of 
these properties for the whole series of alkali hydrides have been demonstrated. Finally, we discussed 
the importance of the accuracy of the representation of the hydrogen electronic affinity in the long- 
range dynamics of the molecular excited states. In particular, we suggested that the empirical model 
based on a diabatization procedure by Gadea and coworkers may not address this issue in the same way 
as our calculation using an extended basis set for the hydrogen atom. Recent developments may soon 
bring comparison between theory and experiment into reach. For instance, an intense supersonic beam 
of LiH molecules has been set up in the perspective of a Stark deceleration experiment [100]. Also, 
very precise Stark spectroscopy has been demonstrated in an ensemble of ultracold KRb molecules in 
their ground state, yielding accurate values for permanent and transition dipole moments [101]. 
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Fig. 8. Transition dipole moments (upper panel) and permanent dipole moments (middle panel) for the five lowest 
1 E + states of RbH computed in the present work (full lines), and by Khelifi et al. [29] without (with) empirical 
correction for hydrogen electron affinity (dashed and dotted lines, respectively). Related potential curves computed 
in the present work are displayed in the lower panel. The link between the presence of avoided crossings and the 
abrupt variations of the dipole moment is highlighted with vertical dashed lines. 
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X 2 E+ 


A 2 E+ 


Ion 


Reference 


Re 


D e (eV) 


Re 


D e (eV) 


LiH+ 


This work 


4.11 


0.1426 


7.44 


0.4994 




Rosmus and Meyer [50] 


4.12 


0.13 








Fuentealba et al. [47] 


4.12 


0.14 








Alikacem and Aubert Frecon [43] 


4.08 


0.13 


7.35 


0.50 




Vojtikrfa/. [51] 


4.15 


0.13 


7.27 


0.41 




Berriche et al. [49] 


4.11 


0.141 


7.45 


0.497 




Dalgarno et al. [1] 


4.13 


0.1404 


7.39 


0.48 




Kimura et al. [37] 






7.47 


0.49 




Magnier [40] 


4.15 


0.13 


7.46 


0.495 


NaH+ 


This work 


4.85 


0.083 


7.88 


0.4594 




Valance [45] 


5.8 


0.02 


8.2 


0.34 




Rosmus and Meyer [50] 


4.90 


0.14 








Melius et al. [46] 


5.1 


0.061 








Olson et al. [52] 


5.1 


0.061 


7.98 


0.469 




Liu et al. [53] 


5.1 


0.064 








Fuentealba et al. [47] 


4.86 


0.08 








Watanabe et al. [39] 


4.65 


0.103 


7.8 


0.435 




Kubach and Sidis [36] 






8.7 


0.39 




Kimura et al. [37] 






7.89 


0.46 




Allan [38] 






7.88 


0.46 




Magnier [41] 


4.9 


0.0615 


7.86 


0.4602 


KH + 


This work 


5.66 


0.062 


8.68 


0.6481 




Valance [45] 


10. 


0.025 


8.6 


0.49 




Fuentealba et al. [47] 


5.88 


0.05 








Olson et al. [52] 


5.66 


0.054 


8.74 


0.614 




Melius et al. [46] 


6.5 


0.022 








Kubach and Sidis [36] 






9.4 


0.49 




Kimura et al. [37] 






8.68 


0.68 




Watanabe et al. [39] 


5. 


0.136 


8.45 


0.625 




Magnier [42] 


5.52 


0.0476 


8.58 


0.6607 


RbH+ 


This work 


5.82 


0.6824 


8.58 


0.6824 




Valance [45] 


13.4 


0.012 


8.8 


0.50 




Von Szenptaly et al. [48] 


5.83 


0.06 








Kubach and Sidis [36] 






10.4 


0.55 




Kimura et al. [37] 






8.87 


0.73 


CsH + 


This work 


5.98 


0.0809 


9.36 


0.793 




Valance [45] 


17.6 


0.03 


10.4 


0.68 




Von Szenptaly et al. [48] 


6.14 


0.03 








Scheidt et al. [35] 






10±1 


0.77± 0.05 




Kubach and Sidis [36] 






10.4 


0.71 




Kimura et al. [37] 






9.20 


0.86 



Table 1. Equilibrium distances R e (in a.u.) and binding energies D e for the X 2 E + and A 2 E + states of alkali 
hydride cations obtained in the present work, compared to previous predictions. 
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This work 








other works 






state 


He 


L>e 


UJe 


ref. 


T) 

He 


D 

L> e 


U)e 


LiH 


X 


3.002 


20166.8 


1398 


[75]+ 


3.015 


20287.7 ±0.3 


1406.9 












[27] 


3.007 


20174 














[74] 


3.0 


20123 


1391 


NaH 


X 


3.54 


15671.3 


1163 


[54]+ 


3.566 


15900±100 


1171.5 












[74] 


3.254 


16050 


1163 


KH 


X 


4.168 


14500.0 


952 


[54] + 


4.23 


14772.7±0.6 


986.6 












[28]" 


4.19 


14365 














L28J 


a 1 n 


1 AHZ(\ A 

i4/3(J.4 














[74] 


4.2 


14937 


961 


RbH 


X 


4.395 


14098.8 


885 


[54] + 


4.47 


14580±600 


937.1 












[29]" 


4.40 


13940 














[29] b 


4.40 


14323.2 














[74] 


4.367 


14800 


912 


CsH 


X 


4.464 


15019.3 


802 


[54] + 


4.71 


14791.2±2 


891.2 












[30]" 


4.47 


14878 














[30f 


4.48 


15253.4 














[74] 


4.626 


14897 


881 



Table 2. Equilibrium distances R e (in a.u.), harmonic constants iu e and potential well depths D e (in cm" ) of the 
ground state of alkali hydrides, compared to available experimental values, and recent theoretical determinations. 
The (a) and (b) labels refer to the calculation without, and with the correction for hydrogen electronic affinity, 
respectively (see text). Note that ref. [54] is a review paper which collected results from many experiments, which 
explains the apparently huge variations of the error limits. 
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This work 








other works 






state 


Re 


D e 




ref. 


Re 






LiH 


A 


4.82 


8698 


241 


[27] 


4.87 


8689 














[55] 


4.91 


8685.6 ±0.3 


244 [75] 




C 1 


3.83 


1267 


390 


[27] 


3.825 


1276.8 






c 2 


10.15 


8361 


390 


[27] 


10.21 


8444 






E 


5.31 


3727 


516 


[27] 


5.358 


2839 






(5) 


4.1 


1455 


486 


[27] 


4.094 


1574 






B 


4.52 


251 


243 


[56] 


4.5 


288.9±0.2 




NaH 


A 


5.97 


10045 


318 


[76] 


6.01 


9997 














[63] 


6.186 


9701 


333.1 [78] 












[54]+ 


6.03 


10143 


317.56 




C 1 


4.45 


661 


432 


[76] 


4.48 


612 






c 2 


11.75 


6501 


219 












E 


5.23 


3709 


683 


[76] 


5.24 


4738 






(5) 


4.57 


1121 


432 


[76] 


4.63 


3927 






B 


5.1 


282 


214 










KH 


A 


6.95 


8884 


243 


[29] 


7.05 


8946.3 














[79] 


7.18 


8710 


288 












[80] 


6.944 


8590 


245.9 












[57]+ 


7.11 


8698 


222.74±0.16 




C 


13.3 


6471 


163 


[29] 


13.63 


6584.5 






E 


5.59 


970. 


232 


[29] 


5.65 


873.5. 






(5) 


5.27 


1028 


362 


[29] 


5.28 


1026.9 






B 


5.84 


294 


186 


[29] 


5.34 


801 




RbH 


A 


7.17 


9053 


217 


[28] 


7.3 


8711.2 














[80] 


6.92 


8710 


232 












[54]+ 


6.270 


8941 


211.74 




C 1 


5.61 


3568 


215 


[28] 


5.65 


3545.9 






c 2 


14.1 


5387 


197 


[28] 


14.4 


5509 






E 


5.74 


553 


245 


[28] 


5.77 


504.3 






(5) 


5.5 


1162 


353 


[28] 


5.6 


1133.3 






B 


5.89 


355 


213 


[28] 


5.92 


377.3 




CsH 


A 


7.31 


8825 


192 


[30] 


7.55 


8851.9 














[66] 


7.52 


7500 














[81] 


8.10 


8630 


196 












[67] 


7.46 


7767 


196 












[54]+ 


7.52 


8130 


196 




C 1 


5.63 


2543 


301 


[30] 


5.63 


2462.5 






C 2 


14.8 


2175 


111 


[30] 


29 


3785 














[81] 


6. 


2180 


242 












[67] 




1930 


144 




E 


5.83 


659 


264 


[30] 


5.86 


642.4 






(5) 


5.69 


1288 


332 


[30] 


5.72 


1354.36 






B 


6.01 


576 


228 


[30] 


6.04 


551.89 





Table 3. Main spectroscopic constants of excited singlet electronic states of alkali hydrides, for the four lowest 
states of 1 E + symmetry (A, C, E and (5)) and the lowest state (B) of 1 n symmetry. The exponents 1 and 2 
labeling the C state refers to the inner and the outer potential well, respectively. Equilibrium distances Re are in 
bohr radii, depth of potential wells D e and harmonic constants u> e in cm -1 . Comparison is provided with available 
experimental results, and with a selection of recently published calculations (mainly from the series of papers of 
Gadea and coworkers. 



©2009 NRC Canada 



Aymar, Deiglmayr, and Dulieu 



23 







a zz 


OCxx 


a 


7 


Hz 


LiH 


This work, for v — 


26.8 


30.3 


29 A 


-3.5 


5.87 




This work 


26.0 


29.8 


28.6 


-3.8 


5.82 




Gready et al. [88] 


34.1 


32.6 


33.1 


1.5 


5.90 




Lazzretti and Zanasi [86] 


21.93 


25.3 


24.2 


-3.4 


6.0 




Roos and Sadlej [62] 


26.3 


29.3 


28.3 


-3. 


5.90 




Bishop and Lam [90] 


26.3 


29.7 


28.6 


-3.4 


6.0 




Vrbikefa/. [92] 


24.6 


30.9 


28.8 


-6.3 


5.76 




Maroulis and Bishop [89] 


21.8 


25.1 


24. 


-3.3 


6.0. 




Sasagane et al. [87] 


22.9 


29.4 


27.3 


-6.4 






Rerat et al. [91] 


26.9 


30.8 


29.5 


-3.9 


5.87 




Sadlej and Urban (a) [93] 


26.7 


30. 


29 


-3.2 


5.909 




Sadlej and Urban (b) [93] 


26.81 


30.01 


28.94 


-3.2 


5.906 




Wharton et al. [96] (exp.) 










5.882 ±0.003 


NaH 


This work, for v — 


55.3 


39.2 


44.5 


16.1 


6.41 




This work 


53.7 


38.8 


43.7 


15 


6.39 




Sadlej and Urban (a) [93] 


54.1 


39.9 


44.6 


14.3 


6.67 




Sadlej and Urban (b) [93] 


58.9 


39.7 


46.1 


19.2 


6.38 




Diercksen and Sadlej [94] 


50.8 


37.7 


42.0 


12.5 


6.44 




Dagdigian [97] (exp.) 










6.4±0.07 


KH 


This work, for v = 


72.0 


44.6 


53.7 


27.4 


8.11 




This work 


71.8 


44.5 


53.6 


27.3 


8.14 




Sadlej and Urban (a) [93] 


66.0 


48.4 


54.3 


17.6 


8.48 




Sadlej and Urban (b) [93] 


85.0 


50.5 


62.0 


34.5 


8.08 


RbH 


This work, for v — 


85.6 


44.1 


57.9 


41.5 


8.36 




This work 


85. 


44.2 


58.0 


41.3 


8.42 




Sadlej and Urban(a) [93] 


71.4 


48.2 


56.0 


23.1 


9.04 




Sadlej and Urban (b) [93] 


95.34 


51.5 


66.1 


43.8 


8.41 


CsH 


This work, for v = 


91.6 


34.7 


53.7 


56.9 


8.22 




This work 


92.6 


37.5 


55.9 


55. 


8.67 



Table 4. Polarizabilities (in a.u.) and permanent dipole moments (in Debye) for the ground state of alkali hydrides, 
either computed for the v — level, or taken at the experimental equilibrium distance, compared to published 
theoretical and experimental data. Values of ref. [88] corresponds to their "CI" calculation, in ref. [87] to their 
"MCTDHF[DQJ" calculations, in [94] to their "MBPT(4)" calculations, and in ref. [93] to their "T(CCSD)" (a) in 
the valence-shell approximation and (b) in a model including correlations from both valence and next-to valence 
shells. 
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